Subclasses in Mixing Correlation-Growth Processes with Randomness 



O 

o 

< 



A. Kolakowska and M. A. Novotny 
Department of Physics and Astronomy, and Center for Computational Sciences, 
P.O. Box 5167, Mississippi State, MS 39762-5167 
(Dated: February 2, 2008) 

We show that in the construction of continuum equations for competitive growth processes that are 
a mixture of random deposition and a correlation process, a distinction must be made within a single 
universality class between depositions that do and do not create voids in the bulk. Within these 
subclasses the bulk morphology is reflected in the surface roughening via nonuniversal prefactors in 
the universal scaling of the surface width. 
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Many complex systems are studied by mapping onto 
a suitable nonequilibrium surface-growth problem 0, |3| ■ 
The dynamics of the buildup of the correlations in a sys- 
tem can be then explored with surface-growth method- 
ologies. Large-scale properties are described within a 
continuum model by universal stochastic growth equa- 
tions and tested with simulation models. The trouble is, 
simple atomistic models are often not adequate to repro- 
duce the complex physics of the observed surface phe- 
nomena that may involve contributions from several uni- 
versal processes, but the continuum description of such 
multicomponent growth has not yet been developed. A 
representative example comes from an applied model in 
computer science |3l when the asynchronous dynamics 
of conservative updates in a system of parallel proces- 
sors is modeled as a virtual-time surface that represents 
nonequilibrium processes in this system. When the load 
per processor is minimal, this dynamics is Kardar-Parisi- 
Zhang (KPZ) universal process [3|. But when the load 
is increased to reflect real operations, the realistic dy- 
namics is a competitive growth that combines a univer- 
sal KPZ process with a random deposition (RD) pro- 
cess, i.e., is of the type RD-f-KPZ 4]. Consequently, it 
is important to know how the nonuniversal properties of 
the component-processes affect the universal scaling of 
a RD-f KPZ process. This is a still unsolved problem of 
nonequilibrium surface growth science. In applied model- 
ing, even if not explicitly assumed, competitive dynamics 
naturally arises. Studying these systems should also con- 
tribute to the understanding of differences between the 
expected and the actual scaling of rough interfaces, often 
encountered both in simulations and in experiments. 

In this work we investigate a connection between sur- 
face roughening and the bulk morphology formed dur- 
ing the deposition in competitive growth RD-f X, where 
X is a correlation growth process of universal dynam- 
ics different from RD. This connection has been already 
established in simulations of competitive growth mod- 
els 13 and binary growth of thin films 0. A new as- 
pect of our study is a direct theoretical link between 
nonuniversal properties of process X, as read from the 
bulk, and the continuum equation that reflects the ob- 



served scaling laws for the competitive RD-fX processes. 
We derive from first principles a continuum equation to 
show that model dependent coefficients do reflect the 
bulk structure. This will lead to a distinction between 
void-producing and desorption growth processes and sim- 
ple absorption processes. As discussed later, this division 
into subclasses is a necessary first step towards a theory 
of many-component processes. In particular, it explains 
variations in scale-dilatation observed in RD-I-X models 
0, 0, 0, S Oi Q • In our analysis we use as an ex- 
ample universal RD-I-KPZ growth processes in (1 + 1) 
dimensions, and generalize our approach to other pro- 
cesses in (1 + n) dimensions. 

A scaling hypothesis for competitive RD-I-X processes 
states that if a correlation growth X occurs with a 
constant probability p, its continuum equation must be 
invariant under the scaling 



h/g{p) , t t/f{p) , 



(1) 



where g{p) and f{p) are arbitrary suitable functions of 
p 6 (0; 1]. This invariance implies that f{p) — cj^{p), and 
when X=KPZ it leads to the KPZ equation J_3J for the 
RD-fKPZ mix 0: 

ht = yof{p)h^. + {\omf'\p)hl + vix, t) , (2) 

where h = h{x, t) is the height field; x and t are the 
spacial and time coordinates, respectively; subscripts de- 
note partial derivatives; r]{x,t) is the white noise; and, 
i/Q and Ao are constants. When Aq = 0, Eq. (|2Jl is the 
Edwards- Wilkinson (EW) equation Q when X=EW. 
When vq = Xq — 0, Eq. JSJ defines universal RD dynam- 
ics. Many simulation models of RD-I-EW and RD-f-KPZ 
growth processes 0, ^3 suggest g{p) — p^ in Eq. 
which leads to the Family- Vicsek universal scaling [1 
of the average surface width t) 0: 



(3) 



For substrates of size L, F(y) describes two limit-regimes 



of evolution: F{y) 



Oijz 



if y ^ 1 (growth); and. 



F(y) ~ const if y ^ 1 (saturation). In Eq. a and 
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z are the universal roughness and dynamic exponents, 
respectively, of the universality class of the correlation 
growth X. The scale-dilatation exponent 8 in scaling pref- 
actors in Eq.Q, however, is nonuniversal. It has been 
observed that in some models 5 sa 1 across universality 
classes, and in some other models < S ^ I within a 
single universahty class Also, there are models 

where prefactors in Eq.lPJ do not at all obey a power law 
m p M- Here, we shall establish that this variation is 
not accidental, but rather reflects the properties of the 
bulk of the deposited material. 

Consider aggregations where particles fall onto a sub- 
strate of L sites, where they may be accepted in accor- 
dance to a rule that generates correlations among the 
sites. This correlation growth occurs with probability p 
and competes with RD growth that occurs with proba- 
bility q = 1—p. When a particle is accepted at a site, the 
site increases its height by Ah. If, e.g., component 1 is 
RD, and component ^ is a correlation growth in the KPZ 
universality class, their corresponding growth equations 
are 

hi,t = m{x,t) , (4) 

h2,t = Voh2,xx + {>^o/'2)hl^^+T]2{x,t), (5) 

where /i„(x,t), n — 1,2, is the column height at x after 
time t when the component n acts alone. Assume for sim- 
plicity that the noise terms are of the same strength, i.e., 
•q = rji — ri2- In two-component growth, when both com- 
ponents act simultaneously together, the column height 
h{x,t) is incremented due to either of the components 
with their corresponding probabilities p and q, p + q — 1: 

Ah{x,t) ^pAh2{x,t) +qAhi{x,t). (6) 

Here, probability p (or q) is the fraction of contributions 
to h from component 2 (or 1). For some processes this 
fraction is identical to a fraction of times when h{x) is in- 
cremented due to component 2 (or 1) for the times from 
to t. However, as explained later, this is not so for 
all processes. In Eq. (jS)), Ahn is understood as "being 
incremented due to the process n," n = 1,2. In this 
statistical sense, Eq. (0 expresses a simplectic decompo- 
sition of Ah{x,t) into its vertex-components A/i„(x,t). 
Dividing Eq. 10 by At, and taking the limit At — > 0, 
gives the equation for time rates, ht — p^2,t + qhi^t^ to 
which we substitute Eqs. Q-©: 

ht = voph2,xx + (Ao/2)p/i2,x + iP + t) . (7) 

In Eq.lO, h{x,t) is the column height that rises at x as 
the result of two processes acting simultaneously from 
the beginning to time t. Here, h2{x,t) is the part of 
h{x, t) that was created by the component 2 in this time. 
The other part was created by component 1. In other 
words, h2{x,t) is so far an unknown fraction of h{x,t). 



To find a relation between h and /i2, one must consider 
nonuniversal properties of aggregation processes. 

We distinguish between the following two groups of 
surface growths. In one group we place all simple ab- 
sorption processes with conserved flux that do not cre- 
ate voids in the bulk of the deposited material. We call 
this group absorption-bulk-compact (ABC) growths. The 
other group, which we call dense- or-lace-bulk (DOLE) 
growths, contains processes that are not ABC-type. The 
DOLB group includes desorption processes that may lead 
to a dense bulk as well as absorptions that lead to the 
formation of voids. Note, RD processes are ABC growth 
processes. As we show in the next paragraph, when com- 
ponent 2 is of the ABC-type, p and q in Eq. ((HJ express 
fractional contributions to h in terms of times, and then 
h2{x, t) = ph{x, t). This is not true when component 2 is 
a DOLB growth. 

Consider a discrete representation of events at coor- 
dinate X. Suppose, there are t deposition events in to- 
tal, with ti events due to component 1, and t2 events 
due to component 2, t ^ ti + t2- In ABC growth, af- 
ter t events, the total column height is h = tAh, where 
contributions from components 1 and 2 are, respectively, 
hi = tiAh and /12 ~ t2Ah. Thus, hi/h — ti/t — q and 
^2/^ = ^2/^ — P- Therefore, in ABC growth h2 — ph, 
and in Eqs. ®-0 P = P and q = q. 

Next consider that the component ^ is a DOLB growth 
that creates voids. Now, an individual deposition event 
due to component 2 not only increases h by A/i, but 
may also result in the creation of voids. The net effect 
is as though component 2 deposited Ah and the voids. 
Therefore, in t2 events, its contribution to the column 
height is h2 — (^2 + m)Ah, where mAh reflects the in- 
crease in height due to the presence of voids. The com- 
ponent 1 is RD, i.e., ABC-type, and hi = tiAh. Af- 
ter t events, the net column height is h — hi + /i2 — 
(t -I- ■m)Ah. Thus, hi/h ti/{t + m) < ti/t = q 
and /12/ft- = {t2 -\- m) / {t -\- m) > t2/t — p. Fractions 
Qoff = hi/h and pcs = h2/h are the effective probabili- 
ties of deposition events due to components 1 and 2, re- 
spectively, as they would result from measuring the col- 
umn height. For some types of two-component growth 
with RD, the probability PcS can be expressed approxi- 
mately as Peff = P^ 01 • F'o'^ DOLB growths with voids 
5 < 1 (because PeS > p)- When the component ^ is a 
DOLB growth with desorption, in the above reasoning 
one should change to ^ —to. This will give geff > q 
and poff < P, and Poff = P^ with 5 > 1. The value of 5 
is specific to the particular deposition process of compo- 
nent 2. Therefore, in DOLB growth /i2 = Pcsh, and in 
Eqs. ®-lIZl) p — PcS and q — qcB- 

In general, h2{x,t) = pesh{x,t) and p = PeS, where 
PcS = p if the correlation component is an ABC growth. 
When the correlation component is a DOLB growth, and 
when the effective probability is well approximated by a 
power law p^ , the result can be summarized as PcS = P^ , 
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where S = 1 for ABC growths and S 1 for DOLB 
growths. This result is combined with Eq. to give the 
continuum equation for the RD+KPZ mix: 



+ {\o/2)p^'hl+7i{x,t) 



(8) 



When in Eq. ISJ Aq = 0, the analogous reasoning gives 
the RD+EW dynamic: 



(9) 



Both results, Eqs.(|SJl-l|njl, are in accord with our for- 
mer derivation that lead to Eq. (01 7] . Matching Eq. ^ 
with Eq. 121 gives f{p) = p^^ , which form of f{p) was 
used formerly to derive the approximate prefactors in 
Eq.©. The inverse of the scaling Q when applied to 
Eqs.®-® transforms them to continuum equations for 
a "pure" correlation processes oi p = 1. Explicitly, it 
collapses all evolution curves w{p,t) (for all L and p) ei- 
ther onto w{l,t) or onto a neighborhood of w{l,t) 0, 
following Eq.®, provided the effective probabilities PcS 
can be approximated by the power-law p^ . When such a 
fit is not possible Eq.Q is obeyed but then the scaling 
prefactors must be expressed directly in terms of effective 
probabilities. This is because the factor p^ in the coeffi- 
cients of Eqs.(|SJ)-l|njl is only a fit to the effective probabil- 
ity PcS- In fig^we give an example of the exact scaling 
where nonuniversal prefactors in Eq.Q are directly ex- 
pressed by pctf via the substitution — > pfff (p) — \/ f{p) 
for the RD+BD model when BD is the NN sticking rule 
[l7l |. Here, the effective probability depends on both p 
and the mean compactness c{p) of the bulk formed in 
the RD+BD process: PcS = 1 — qc(j>) The perfect 

data collapse in the full range of p G (0; 1], seen in fig^ 
can be contrasted with fig. 5 of Ref.|3 that shows only 
an approximate data collapse for the same system with 
the best fit exponent S « 0.41 in Eq.©. It needs to be 
said explicitly that the scaling where 5 = 1/2 in Eq.©, 
proposed in Refs.0, 0] for RD+BD models, docs not 
produce data collapse at all. The RD+BD model when 
BD is the NNN sticking rule provides an example 
where Pcs{p), and thus the nonuniversal prefactors f{p) 
and g{p) in Family- Vicsek universal scaling, cannot be 
expressed by a power-law p^ . In this system the sur- 
face roughening obeys power laws in effective probability 
that incorporates either the compactness or the voidness 
of the bulk, which gives excellent data collapse of w{p, t), 
similar to that seen in figH |ld| • 

The approach introduced here by the example of KPZ 
processes, can be applied to a broad range of stochastic 
growth models RD+X, where component 2 can be any 
isotropic growth in (1 + n) dimensions: 



h2.tix, t) = GQ12) + mix, t) 



(10) 



where x is n dimensional, and the operator G represents 
only local interactions Il7l|. In the general case, Eq.© 
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FIG. 1: (color on line) Scaled time-evolution w^{p,t) in the 
RD+BD model. In this example, the scaling function f{p) = 
(j>) explicitly incorporates the compactness c(p) of the bulk 
formed in the RD+BD process. Here, L = 500, 2a — 1, and 
averaging was performed over 400 configurations. 



is written as hf = Poff^2,t + QcS^i.t, and combined with 
Eqs.Q) and ifTUI) . to find for the competitive growth 



ht{x, t) = pcsG{pc«h) + ri{x, t) . 



(11) 



where 77 = (1 — Pcs)r\\ + PcbV2, and the noise strengths 
may be different. Eas. H10|l - Hll|) represent the same uni- 
versality class since the multiplication by poff does not 
modify local interactions: PcS affects the noise strength 
and the gradient of the height field, but it does not gener- 
ate new terms other than those already given by operator 
G. Hence, if a correlation growth belongs to a given uni- 
versality class, its mix with RD will remain in the same 
class. Elementary calculations show that Ea. lll|l is in- 
variant under the scaling g{p)h{x,t) = h'{x,t' = f{p)t). 
If gip) = Pcsip) and /(p) = p^gip), and if the noise 
strengths are the same, this scaling maps the univer- 
sal dynamics Hll|) of RD+X onto the universal dynamics 
of X. In this case the invariance implies g(j))w{p,t) = 
w'{f{p)t), where w'[-) has universal scaling properties of 
the process X. When X is either in the KPZ or in the 
EW universality class, and if additionally pcff ~ P*, we 
recover Eq.©. 

When both the RD and the correlation component 2 
have deposits of unit height, when peff ~ p'^, we have 
(5 = 1 if component 2 is of the ABC-type; and, (5 ^ 1 if 
it is of the DOLB-type. In the latter case, the value 
of the exponent 5 is specific to component 2. When 
Pcff incorporates explicitly bulk properties, the scaling is 
g{p) = p^fj(p), where the new scale-dilatation exponent i5 
is obtained from the slope of hi-up'ip) plotted vs Inpeff(p) 
at saturation. In DOLB growth with voids, peff can be 
determined by measuring the mean density of voids in 
the bulk (Fig^l. Similarly, in DOLB growth with des- 
orption, poft is connected to the mean fraction of the 
removed material |l6l| . 

The analysis presented here explains scaling results of 
the followin g rn ixed-growth models in (1 + 1) dimensions. 
Model A I3, lid lUl : component 2 is RD with surface re- 
laxation. Model B 0: component 2 simulates a deposi- 
tion of a sticky non-granular material of variable droplet 
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FIG. 2: Scaled widths at saturation w vs the parameter ■ 
(a) and (b) are for Models A and B, respectively; (c) and (d) 
are for Models C and D, respectively. Reference lines have 
slope 1. Data are scaled with the a values shown. 



size. Model C 0, [T3, El : component 2 is the NN sticking 
rule of BD. Model D |3,[3- component ^ is a deposition 
of Poisson-random numbers to the local surface minima. 
Models A and B are ABC growths in the EW universal- 
ity class, where poff — P and 6 — 1 (Fig|2J). Models C 
and D belong to the KPZ universality class. Model C is 
an example of DOLB growth with voids, with a 53.2% 
void density in the bulk when p — I, and in this case 
S « 0.41 < 1. Model D is a DOLB-type growth that 
produces a compact bulk but component 2 is flux non- 
conserving, and here 6 ~ 1. Extensions of Models A 
and C to {1 + n) dimensions |U, n = 2,3, yield results 
that conform to our theoretical predictions of poft ~ 
with (5 7^ 1 for mixing RD with DOLB processes, and 
S — 1 for mixing RD with ABC processes. Additional 
examples include cases when component ^ is a restricted 
Kim-Kosterlitz solid-on-solid model [s^ , and when it rep- 
resents the Villain-Lai-Das Sarma universality class [l9j . 

The extension of the approach presented here to other 
competitive growth processes may provide a tool to un- 
derstand the observed dynamics of surface growth. Re- 
alistic systems may involve many component-processes, 
some of which may be dominant. Within our formalism 
a departure point may be a generalization of Eq. (|5J : 



(12) 



where the summation is over contributing processes, and 
Ah^'''> is the column-height increment due to the fcth pro- 
cess. In first approximation component-processes are not 
explicitly correlated. Each process is encountered with 
probability pfc, X^fcPfc ^ 1' and contributes to the growth 
with an effective probability p^^^ X^fcPiff ~ 1- the 
trivial case of all components being ABC type models 



(k) 

with unit mean deposit height = pk- For a DOLB 

(k) 

growth pjjjj will have to be determined. Depending on the 
model, this can be done by analyzing the growth when 
process k acts alone, and measuring either the mean bulk 
density or the mean fraction of the detached material or 
both 16]. Simplectic decompositions like the one pro- 
posed in Ea. H12l) have a long history of applications in 
many diverse fields. 

In summary, the derived continuum equations and the 
resulting scaling show that model-dependent prefactors 
in universal scaling laws can be determined from bulk 
structures. This necessitates the distinction between 
the absorption-bulk-compact and the dense-or-lace-bulk 
growth processes in the analysis of competitive mixed- 
growth models. 
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